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Abstract 

Theoretically there are two approaches to predict single spin azimuthal asymmetries. One is to 
take transverse momenta of partons into account by using transverse momentum dependent parton 
distributions, while another is to take asymmetries as a twist-3 effect. The nonperturbative effects in 
these approaches are parameterized with different matrix elements and predictions can be different. 
Recently, gauge invariant definitions of transverse momentum dependent parton distributions were 
derived. With these definitions it can be shown that there are relations between nonperturbative 
matrix elements in two approaches. These relations may enable us to unify two approaches and to 
have unique predictions for single spin azimuthal asymmetries. In this letter we derive these relations 
by using time-reversal symmetry and show that even with these relations the single spin azimuthal 
asymmetry in Drell-Yan process is predicted differently in different approaches. 

PACS numbers: I3.85.Qk, I3.88.+e 
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Single spin azimuthal asymmetry provides a new tool to study structure of hadrons because the 
asymmetry is sensitive to correlations between quarks and gluons as partons inside a hadron and to 
orbital angular momenta of these partons. Experimentally, such an asymmetry was observed for inclusive 
production of pion in polarized proton antiproton scattering with center-of-mass energy y/s = 20GeV by 
E704 collaboration 1 J. The asymmetry is large for charged pion, while for 7r° production it is consistent 
with zero when the transverse momentum kx is smaller than 3GeV, and it tends to a positive value when 

becomes large. In semi inclusive deep-inelastic scattering(SIDIS) significant asymmetries were also 
observed in production of pion and kaon by HERMES 2 . Asymmetries in polarized proton scattering 
are currently studied by STAR at RHIC. Large spin effects are observed in their preliminary results after 
first run. In SIDIS measurements of asymmetries with a transversely polarized target were reported by 
the SMC collaboration [2]. Experiments with a transversely polarized target are now under study by 
HERMES and COMPASS @|E|. 

Single spin azimuthal asymmetry is a T-odd effect and helicity-flip amplitudes are involved. Per- 
turbatively T-odd effects can be generated at loop-level in hard scattering of active partons of hadrons. 
Because the quark-gluon coupling of QCD conserves helicities in the massless limit, the T-odd effects are 
proportional to quark masses which can be neglected. Therefore the observed T-odd effects can not be 
explained by those T-odd effects arising from hard scattering and are related to nonperturbative nature 
of hadrons. Indeed, these T-odd effects can be generated from final- or initial interactions between active 
partons involved in the hard scattering and remnant partons in hadrons |15l 1231 12*1] . The effect of these 
interactions can be represented by gauge links in definitions of parton distributions [231 124j . Theoretically 
there are two approaches to explain single spin azimuthal asymmetry by taking nonperturbative nature 
of hadrons into account. One is to take transverse momenta fcy of partons in a hadron into account 
where one uses transverse momentum dependent parton distributions to parameterize nonperturbative 
effects. For a polarized hadron as an initial state the effect is parameterized by Sivers function 6 , while 
for a hadron observed in a final state the T-odd effect related to this hadron is parameterized by Collins 
function^. For semi-inclusive deep inelastic scattering, both functions can make contributions to the 
observed single spin azimuthal asymmetry. Single spin azimuthal asymmetry has been studied in terms 
of these functions [SJ OH E3 CCD EE21 EES] • These functions have been also studied with models [TH IT51 fTTH ITTj . 
Another approach, called Qiu-Sterman mechanism, is that the T-odd effect is produced by taking twist-3 
effect into account and it is proportional to quark-gluon correlations inside a hadron ^Hj. The fact that 
T-odd effects can be generated at twist-3 level was also pointed out in ^2]- This approach was used to 
make predictions for various processes in |181 12U| . It is interesting to note that at first look the physical 
reason for the effect is different in different approaches. In the first approach the helicity of a initial 
hadron is changed because of orbital angular momenta of partons. This can be seen clearly in terms of 
light-cone wave functions |21|. In the second approach the helicity flip is caused by nonzero spin of the 
gluon which is correlated with other partons. Predictions based on different approaches are different. A 
question arises why there are two physical origins for one effect? 

This question has been answered partly by recent studies of transverse momentum dependent parton 
distributions [221 1231 124j . which are involved in the first approach. It has been shown that gauge links in 
these distributions play an important role to incorporate T-odd effects introduced by final state interac- 
tions. In particular, additional gauge links should be included in the definitions of these distributions 24 . 
With these gauge links it is possible to relate the second hr moment of Sivers function to the twist-3 
matrix element in the second approach ,25]. With relations between nonperturbative matrix elements in 
different approaches it may be possible to unify two approaches and to have unique predictions for single 
spin azimuthal asymmetries. In this letter we will show that predictions based on two approaches are 
still different, although such relations exist. We will show this in detail with Drell-Yan process. Before 
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showing this we give another derivation of relations between second moments of T-odd distributions 
and twist-3 matrix elements by using time-reversal symmetry of QCD. 

We consider a proton moving in the z-direction with the momentum P and the transverse spin St • We 
use a light-cone coordinate system and introduce two light-cone vectors: ra M = (0, 1,0,0), P = (1,0,0,0) 
and n ■ I = 1. Neglecting the proton mass we have P^ = (P + , 0,0,0). Taking transverse momenta of 
partons in quark-quark correlation in a proton into account, there are two T-odd parton distribution 
functions appearing in a Drell-Yan process which take effects of initial state interaction into account. 
They can be defined as 26 

/iW*. = \J^^e- ik <{(p,s T \m^v(om\p^T) 

-(s T -» -s T )}, 

hi jDY (x,k 2 T )kir = -J^^e- ik <-(P\m^ +i V(0m\P), (1) 

with £ M = (0, £~, £t) an d = £p ( 7fj,i/n p l' 7 . The momentum k is k^ = (xP + , 0, k T ). The matrix element 
in the last line is spin averaged. The function f^j, DY (x, kj<) is the Sivers function for Drell-Yan processes. 
V(£) is a product of gauge links to make the matrix element gauge invariant, it takes the effect of initial 
state interaction in Drell-Yan process into account. If one can take V(£) as a unit matrix, then one can 
show with time-reversal symmetry that both correlation functions are zero. It is important to note that 
V(£) is not a unit matrix even in the light-cone gauge n-G = 0, additional gauge links must be introduced 
to make the definitions gauge invariant 24j. We will take the light-cone gauge. In this gauge V(£) reads: 

V(0 = F-oc(£t) = Pexp (igj^ dfr • G T (0,r = -oo,£ T )l . (2) 

This gauge link takes effects of initial state interaction into account and it can be derived in a similar 
way as in SIDIS 24]. The difference is that the gauge link is at £ _ = — oo because it is for initial state 
interaction. 

Under parity- and time-reversal transformation, we obtain for the matrix element: 

(P,s T |^(e)7 + ^-oo(eT)V(0)|P,s T ) = <P,-s t |^(07 + ^tMO)|P,-s t ), (3) 

with 

V^t) = Pexp lig [ T d£ T ■ G T (0,r = D°,£r) I ■ ( 4 ) 



o 

Similarly one can define two T-odd parton distribution functions appearing in deep inelastic processes 
which take effects of final state interaction into account. The two functions f^r dis( x ^t) an d hiDis( x i ^t) 
are defined similarly as in Eq.(l), but with the gauge link V(£) is replaced with VJ>o(£t)- These two func- 
tions are related to those in Drell-Yan processes with time- reversal symmetry. With Eq.( |2J) we can 
write: 

f^ DY (x,k 2 T )e^ u k^ T = \J ^-^ e -^(P,s T |^(07 + [^oo(eT)-U 00 (e T )]^(0)|P,s T ) 

-(s T -> -s T )}, 
1 [d£-d 2 £ T 



hi DY (x,k 2 T )k^ = "2 / ^jf e " ,fc '^l^(^ +i ^-»(?t) - ^(?t)] t*(0)|P). (5) 
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It is expected that the function /j^DY and hj~ DY decrease rapidly with increasing kx- Then T-odd 
effects related to them can be estimated at leading order by the second moment of k F of the left hand 
side in Eq.(5): 

Kf(x) = J d 2 k T k^ DY (x, k 2 T )EL^s v T = ~\e aa s Ta J d 2 k T \k T \ 2 f^ DY (x, k 2 T ), 
K^{x) = J d 2 k T k»hi DY {x, k\)k v T = ^{n^l u + n u P - gT) J d 2 k T \k T \ 2 h^ DY (x, (6) 

With Eq.(|5J) these moments can be expressed in term of matrix elements. Taking KJ as an example, we 
have 

Kf(x) = -I|^ e -^ + «-^{(P,s T |^(07 + [^oo(eT)-V^(eT)]^(0)|P,ST) 

_( ST ^ _ St )}| 5t=0 . (7) 
Carrying out the derivatives we have: 

Kf(x) = i|^e-- p+ «"-{ 5 (P,ST|^(rn)7 + [G?(0,oo,0,0)-G?(0,-oo,0,0)]^(0)|P,s T ) 

-(s T -s T )}. (8) 

Now one can show that K a is related to the twist-3 quark gluon correlation Tp(x,x) introduced in 
1181. The correlation function is defined as: 



T F ( Xl ,x 2 )e^nJ a s Tp = - 9 - J ^ e "^(— ^-m*^ 

■{(P, s T \4>( yi n) 7 + G + »(y 2 n)i>(0)\P, s T > - (s T - -s T )}. (9) 
where we include the coupling constant g into the definition. It is straightforward to obtain: 

T F (x,x) = J d 2 k T \k T \ 2 f^ DY (x,k 2 ). (10) 

Similarly we have: 

T H (x,x) = J d 2 k T \k T \ 2 h^ DY (x,k 2 ,), (11) 
where Th is defined with a twist-3 operator: 

T H (x 1 ,x 2 ) = g J ^^ e ~^fe--OP + -m-iP + ( p|^ (yin)c7 +M G + i(y2n)V;(0) |p ) . (12 ) 

The relations in Eq.(10,ll) clearly show that the effect of orbital angular momenta of quarks is closely 
related to that of quark-gluon correlations because of gauge invariance. These relations also show that 
the nonperturbative effects in the two approaches for single spin azimuthal asymmetries are the same. 
However, it should be noted that perturbative coefficients in these two approaches are calculated in 
different ways. In the first approach one uses /cy-factorization, while collinear expansion is used in the 
second approach. If the perturbative coefficients in the two approaches are related in a consistent way so 
that the predicted single spin azimuthal asymmetries are same, then we may have an unique prediction 
for single spin azimuthal asymmetries and the question asked before is fully answered. It is difficult to 
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establish a general relation between the perturbative coefficients, since they are differently calculated in 
different ways and are different in different processes. But we can show that the single spin azimuthal 
asymmetry in Drell-Yan process is differently predicted between the two approaches. 
Now we calculate single spin asymmetry in Drell-Yan process: 

A(P A , s T ) + B(P B ) -» r (Pi) + Z+(P 2 ) + A, (13) 

where the proton A is transversely polarized with the spin vector st and moves in the +z-direction. The 
x-direction is chosen as the direction of St- S = (Pa + Pb) 2 - The hadron B is unpolarized and moves in 
— z-direction. We will calculate the single spin azimuthal asymmetry at leading orders, where the lepton 
pair has a small transverse momentum. We assume that the solid angle £1(8, <jf) of the produced lepton 
in the center-of-mass frame of the produced lepton pair and the invariant mass Q 2 of the lepton pair is 
observed. The single spin asymmetry is defined as: 



( da(S T ) _ da(-S T ) \ i ( da(S T ) da(-S T ) \ 

V dQ 2 dQ dQ 2 dQ J ' \ dQ 2 dfl dQ 2 dQ )' 1 ' 

The asymmetry is calculated in .27 with Qiu-Sterman mechanism. The result reads: 

A _ 1 sin26»sin0 fJ2q e l I dx A dx B 5(Q 2 - x A x B S)T Fiq / A (xA, x A )fq/ B (X B ) + \ 
Q l + cos 2 6> ^ Y, q e 2 J dx A dx B f q/A (x A )f q -/ B (x B ) ' J 

where T Fq / P (x, x) is defined in Eq.(9). The subscriber q/A denotes the distribution of q in hadron A. It 
should be noted that the summation YJ and also in the below is over all quark and antiquark flavors, 
i.e., q can be an antiquark in the summation. We only keep the term with T Fq / A . The • • • represents 
another term proportional to T Hq i B which is irrelevant in this letter. 

In order to make comparison of two approaches we study the single spin asymmetries with kx de- 
pendent distributions. At tree level the partonic process is just qq — > l + P . The cross section can be 
written: 

2 1 f d^Pi d 3 P 2 T „ rav 1 



a ~ 6 25 J (2tt) 3 2P° (2vr)32P 2 ° ' Q4 ( 16 ) 

where the leptonic tensor L^ v and hadronic tensor are given by: 

• [(PA,ST\qj(0)qi^i)\PA,s T )(PB\qk(0)qi(^)\PB) +•••]> (17) 
where the • • • denotes power-suppressed terms. The quark q and q carries the momentum 

k A = x a Pa + k A T 

k B = x B P B + k BT (18) 



respectively. The quark density matrix with kx dependence 

1 fd£d^r ik< 



* ij (x,b r ;P,S) = -l "3 e^(P,5|^(0)^(g)|P,5)| g+=0 (19) 
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can be parameterized as|26| 

$(x, k T ; P, S) = ^{hf + /li.cye^Vr A£4 + 31*75^ + h XT io s v T 

+hi s ia^j 5 l^ + h{ DY a^l u } (20) 

where the Sivers function is fif B y The function fxTDY an d h 1DY is defined in Eq.(l). We changed 
the notation of [IS] slightly by replacing 1/M with 1. With this parameterization we have 

= - J2 4 [ dk+d 2 k AT dk B d 2 k BT (2ir) 4 5 4 (Q -k A - k B )(k AT x s T ) • z 

■fiT, q /A(kA)fi, q/ B(kBW u ~ l u n» - l»n v ] +■■■, (21) 

where we keep only terms with /f^^y z denotes the direction of the z-axis. It should be noted that 
the total momentum Q of the lepton pair has nonzero transverse components in general. It depends on 
transverse momenta of incoming partons. It is now straightforward to calculate the asymmetry defined 
in Eq.(14). Since the asymmetry is defined as a distribution of variables in the center-mass frame of 
the lepton pair, we need to specify the frame. We assume that the center-mass frame is obtained from 
laboratory frame by a Lorentz boost only. This is conveniently used in experiment. In the center-mass 
frame the lepton / + and l~ has the momentum £; 2 and k\ respectively. The momentum k\ and &2 read: 

VQ 2 

ki = — - — (l,sin0sin</>, sin cos ^, cos 0), 
VQ 2 

fc 2 = — - — (1, — sin 9 sin <p, — s'mO cos (j), — cos (f>) (22) 
The momentum Pj(i = 1, 2) in the laboratory frame is related to ki(i = 1, 2) by the boost: 



VQ 2 " V 1 Q° 

p ' = k ' + ^ Q+ {^- 1 )^ (23> 

The phase space integration is invariant under the boost. Using Eq.(23) one can express L^ u in term 
of kifa and Q. For transverse momentum independent parton distributions one expects in general that 
they decrease rapidly with increasing transverse momenta. Hence an expansion of the perturbative part 
in transverse momenta is an good approximation. Keeping the first non-zero order in the expansion of 
k A T and ksT-, we obtain the asymmetry: 

1 sin20sin0 ^ 2 f 2 x B - x A 

An = q ' TT^e ^ J dxAdXB6(Q - XaXbS) 2{^ + ^? 

■fq/B{x B )- j d 2 k T \k T \ 2 fi TDYtq/A (x A ,kT)/J2 e2 q / dx A dx B f q / A {x A ) f q j B {x B ) H . 

J q 

1 sin 20 sin 2 f 2 x A - x B 

-Q - TT^Y* J dXAdXBS{Q - XAXbS) 2(^x- a + ^) 2 

■fq/B(xB) ■ T F>q/ p(x A ,x A )/J2 e2 q J dx A dx B f q/A (x A )f q/B (x B ) -\ . (24) 



6 



In the above equation we have assumed that the initial hadrons are in a center-mass frame, i.e., P A = P B . 
In the last step we have used the relation in Eq.(lO). Again, the summation J2 q 1S over all quark and 
antiquark flavors. It is clearly that the asymmetry here is different than that in Eq.(15) because of the 
factor (xa — xb) /%{\/xa + \fx~s) 2 ■ If the factor was 1, then the asymmetry would be the same. Hence, 
the asymmetry obtained by two approaches will have the same angular distribution but the normalization 
is different. Since the factor can be positive or negative, the asymmetry from the two approaches can 
even have different signs. 

It should be noted that the hadronic tensor calculated with the parameterization in Eq.(20) is not 
invariant under electromagnetic gauge transformation. This can be seen by evaluating Q^W^ V with W^ u 
given in Eq.(21). The reason is that the partons involved in the hard scattering have nonzero transverse 
momenta and 7 • I = 7" is contracted with the hard part of qq — ► according to the first two terms in 
Eq.(20). With nonzero transverse momenta the contraction with 7" does not make the qq pair on-shell. 
Hence the U em (l) gauge invariance is not preserved. It is also indicated |28j that the decomposition in 
Eq.(20) may need to be reexamined because the density matrix element also depends on the vector 
due to gauge links. This dependence is neglected in Eq.(20). There are possibly many ways to restore the 
gauge invariance. In this letter we simply make the initial parton on-shell by replacing 7" with j-k/n-k 
for the first two terms in Eq.(20), i.e., 

$(s, k T ; P, s T ) = ~{/i(x, k 2 T ) £± + f^ DY (x, k 2 T )l^e ±pa k p T s$ + ■■■. (25) 
4 k ■ n ' k ■ n 

The hadronic tensor obtained with Eq.(25) is ?7 ern (l)-gauge invariant. It is: 

= i e l I dk+d 2 k A Tdk B d 2 k BT (2ir) 4 5\Q - k A - k B )(k AT x s T ) • z 

6 q J 

{kA .n)(k B .l) [rkA " ^ ~ ^ ~ k A k B]flT,DY, q /A(kA)h,,/B(kB). (26) 



It is straightforward to show Q^W^ ~ k 2 A k u B + k 2 B k u A ~ k\. Hence the gauge invariance is preserved up 
to order k\. The asymmetry calculated with this tensor will be gauge invariant. The result of An can 
be obtained from Eq.(24) by replacing the factor (x B — xa) /2(\/x~a + ^/xs) 2 with ^fxsj^^fxA + ^/x~b). 
Therefore, even after we make the hadronic tensor gauge invariant, the obtained asymmetry An is still 
different than that in Eq.(15) from the second approach. It is interesting to see how the same asymmetry 
in Eq.(15) can be obtained by starting from Eq.(26). If we replace the tensor [g^kA • k B — k A k B — k A k B ] 
with [g^k A -k' B -k'Xk%-k A k^} where k% = (0, k B , — kAT?) and neglect the dependence of lepton momenta 
on transverse momenta of incoming partons, we indeed obtain the asymmetry in Eq.(15) with the same 
normalization, but with an extra negative sign. However, the transverse momentum of the lepton pair 
can not be neglected and we can not do the replacement. 

To summarize: There are two different approaches for single spin azimuthal asymmetries. Using time- 
reversal symmetry, we give in this letter an detailed derivation of the relations between fey dependent T- 
odd distributions and twist-3 quark-gluon correlators, which are used in different approaches, respectively. 
These relations show that the physical origin in the two different approaches for single spin azimuthal 
asymmetries is the same because of the gauge invariance. With these relations it may be expected to 
unify these two approaches and to delivery an unique prediction for single spin azimuthal asymmetries. 
We have studied in detail the single spin azimuthal asymmetry in Drell-Yan process with Sivers function 
and found that predictions from different approaches are different even with these relations. The k? 
factorization used for single spin azimuthal asymmetries does not respect the Z7 em (l)-gauge invariance. 
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This problem may be solved by changing the projection of the perturbative part slightly. But even after 
this changing the predicted asymmetry is still different. Our study shows clearly that different approaches 
give different predictions for the same physical effect in Drell-Yan process and one can expect that the 
same situation will also appear in other processes. Therefore, at present we have not an unique prediction 
for single spin azimuthal asymmetries in Drell-Yan process at least and this problem needs to be studied 
further. 
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